一种从经纬度网格资料获得正方形网格资料的插值方案     

邓爱军  卢敬华 《成都信息工程学院学报》1990,(Z1)

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。  相似文献   

Numerical calculations in the orbital determination of an artificial satellite for a long arc     

Lin Liu  Xinhao Liao 《Celestial Mechanics and Dynamical Astronomy》1994,59(3):221-235

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.  相似文献   

A numerical method for differential equations of second order     

Douxing Cui  Tian-Yi Huang 《Celestial Mechanics and Dynamical Astronomy》1993,55(4):405-409

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.  相似文献   

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target     

K. Wünnemann  B.A. Ivanov 《Planetary and Space Science》2003,51(13):831-845

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.  相似文献   

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires     

Stefan Wessling  Claudia Kuenzer  Winfried Kessels  Manfred W. Wuttke 《International Journal of Coal Geology》2008,74(3-4):175-184

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.  相似文献   

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm     

S. Amir Reza Beyabanaki  Roozbeh Geraili Mikola  Kianoosh Hatami   《Computers and Geotechnics》2008,35(3):346-356

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.  相似文献   

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir   总被引：1，自引：0，他引：1  

Karsten Pruess 《Environmental Geology》2008,54(8):1677-1686

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.  相似文献   

Determination of historical emission of heavy metals into the atmosphere: Celje case study     

Gorazd Žibret 《Environmental Geology》2008,56(1):189-196

The research shows that in the Celje area (Slovenia), the historical anthropogenical emissions are 1,712 tons of Zn and 9.1 tons of Cd. For Zn, this value represents approximately 0.3% of the total Zn production in that area. Close to the former zinc smelting plant, the “Zn precipitation” has been estimated to be up to 0.036 mm. The 100-year Zn production left behind a heavily contaminated area with maximum concentrations of Zn of up to 5.6% in attic dust and 0.85% in the soil, and 456 mg/kg of Cd in attic dust and 59.1 mg/kg in the soil. The calculation of historical emissions is based on the data of heavy metals concentration in the attic dust at 98 sampling points and on the data from 19 measurement sites of the weight of total monthly air deposit. The main idea behind determining past emissions is that when the weight of the deposited dust on a small area is multiplied by the concentration of the element in that area, the mass of the polluter which has been transported to the place of interest by air can be calculated. If we sum up all the weight over the whole geochemical anomaly, we get the quantity of historical emissions.  相似文献   

An operational multiscale system for hazards prediction, mapping, and response     

David P. Bacon  Nash’at N. Ahmad  Thomas J. Dunn  Michael C. Monteith  Ananthakrishna Sarma 《Natural Hazards》2008,44(3):317-327

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.  相似文献   

Numerical modeling of tide-surge interaction along Orissa coast of India     

Prakash C. Sinha  Indu Jain  Neetu Bhardwaj  Ambarukhana D. Rao  Shishir K. Dube 《Natural Hazards》2008,45(3):413-427

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.  相似文献